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Pasadena,  Cadifornia,  on  Air  Force  Contract  AF33(616)-7341  under  Project  Nr 
1370,  "Dynamic  Problems  in  Flight  Vehicles",  Task  Nr  13456,  "Resonant  Fatigue 
of  Aircraft  Structures".  The  work  was  administered  under  the  direction  of  Flight 
Dynamics  Laboratory,  Aeronautical  Systems  Division.  Mr.  Maurice  J.  Cote  was 
task  engineer  for  the  laboratory. 

The  studies  presented  began  in  May  I960,  and  were  concluded  in  March  1961. 
The  nominal  authors  performed  the  bulk  of  the  work,  and  other  members  of  the 
NESCO  Technical  Staff  contributed  in  their  specialties  as  required.  H.  C. 
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The  cooperation  of  Boeing  Airplane  Company,  Douglas  Aircraft  Company, 
General  Dynamics  Corporation,  and  Lockheed  Aircraft  Company  is  gratefully 
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ABSTRACT 


The  stress  in  aircraft  structure  resulting  from  multi-mode  response  to 
sonic  excitation  is  resolved  into  an  alternating  stress  component  superimposed 
upon  a  slowly  varying  mean  stress  component.  It  is  then  found  tlial:  the  proba¬ 
bility  distribution  of  those  components  is  nearly  independent  of  the  number  of 
modes  participating  in  the  response.  This  finding  could  have  considerable  ap¬ 
plication  in  simplifying  fatigue  analysis  and  testing  if  it  may  be  shown  that  the 
mean  stress  component  has  only  low  order  effects  upon  fatigue  life.  Some  of 
the  possible  applications  are  presented, 

A  test  progran'  for  establishing  the  significance  of  the  mean  stress  compo¬ 
nent  is  outlined. 
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1.0  INTRODUCTION 


This  report  presents  tlie  results  of  a  series  of  analytical  and  experimental 
investigations  performed  to  determine  the  probability  distribution  of  the  stress 
response  of  structures  exposed  to  high  intensity  acoustic  noise.  The  object  of  the 
investigations  was  to  reduce  the  number  of  parameters  required  to  describe  the 
response  of  structures  which  may  be  subject  to  fatigue  damage  caused  by  acoustic 
noise.  Any  such  reduction  could  lead  to  simpler  and  more  accurate  predictions 
of  fatigue  life. 

Accurate  prediction  of  fatigue  life  of  structures  loaded  by  random  phenomena 
is  presently  conceded  to  be  beyond  the  state  of  the  art.  The  reasons  for  the  diffi¬ 
culty  are  two-fold: 

1.  Techniques  for  accurate  prediction  of  the  local  stress -time  history 

have  not  been  developed  sufficiently. 

2.  Even  when  the  stress-time  history  is  known  with  sufficient  accuracy 

it  is  seldom  possible  to  make  a  direct  comparison  between  the  known  stress¬ 
time  history  and  one  whose  effects  are  known  through  laboratory  testing. 

Neither  of  those  difficulties  is  unique  to  fatigue.  However,  the  latter  diffi¬ 
culty  in  particular  is  more  marked  in  long  life  fatigue  than  it  is  in  high  stress 
loading.  This  is  true  in  spite  of  the  great  amount  of  laboratory  testing  which  has 
been  performed  in  the  past  for  determining  t)ic  fatigue  characteristics  of  materials 
because  the  groat  majority  of  that  testing  has  been  conducted  using  ideal  loading 
conditions  whicli  are  rarely  duplicated  in  real  structure. 

The  more  widely  recognized  difference  between  typical  laboratory  loading 
procedures  and  typical  conditions  of  actual  loading  has  been  tlic  existence  in  many 
practical  structures  of  a  wide  range  of  stress  levels  applied  in  random  sequence, 
whereas  laboratory  tests  have  principally  been  restricted  to  repeated  applications 
to  failure  of  a  single  stress  level.  Many  theories  liave  been  proposed  (References 
1,  2,  3,  4)  for  correlating  the  result.*:  of  these  laboratory  fatigue  tests  with  random 
level  loading,  and  success  may  be  attained  within  the  forseeable  future. 

However,  success  in  that  correlation  may  not  provide  a  sufficient  tool  for 
predicting  life,  because  those  efforts  to  derive  theories  of  cumulative  damage  neg¬ 
lect  an  additional  difference  between  laboratory  loading  procedures  and  the  com¬ 
mon  conditions  of  life.  That  difference  is  illustrated  in  Figure  1,  where  it  is  in¬ 
dicated  that  laboratory  tests  almost  invariably  consist  of  repeated  cycles  of  alter¬ 
nating  stress  about  a  constant  mean  stress.  This  is  true  whether  the  more  com¬ 
mon,  single  level  test  (Figure  1-a)  or  the  two  or  more  level  tests  performed  to 
gain  knowledge  about  the  stress  interaction  (Figure  1-b)  are  used.  Figure  1-c, 
illustrating  a  typical  aircraft  stress  record,  differs  from  the  laboratory  tests  in 
that  it  does  not  have  a  horizontal  axis  of  stres.5  symmetry.  The  significance  of 
the  lack  of  symmetry  has  not  been  established. 

In  the  extreme,  the  lack  of  symmetry  could  be,  as  Weibull  (Reference  5)  has 
suggested,  sufficiently  influential  that  it  will  be  impossible  to  correlate  laboratory 


Manuscript  released  by  the  authors  14  March  1961  for  publication  as  a  WADD 
Technical  Report. 
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tests  v/lth  experience  unless  the  asymmetry  o£  the  aircraft  stress  is  essentially 
reproduced  in  the  laboratory.  If  that  position  is  found  to  be  correct,  then  the 
efforts  being  expended  by  students  of  stress  Interaction  and  cumulative  damage 
are  of  little  value,  since  every  real  structure  will  have  different  asymmetries 
because  of  different  modal  participation  in  the  total  stress, 

Another  approach  to  the  topic  of  dealing  with  asymmetric  stress  patterns, 
suggested  by  Schjelderup  (Reference  6),  involves  recognizing  that  any  stress,  no 
matter  how  complex,  can  be  transformed  into  a  superposition  of  alternating 
stresses*  and  mean  stresses.*’*'  With  that  transformation,  the  alternating  stresses 
become  symmetrical  by  definition.  At  the  same  time  the  variance  of  the  mean 
stress  frequently  becomes  so  small  relative  to  the  variance  of  the  alternating  stress 
that  it  may  be  approximately  correct  to  replace  the  mean  stress  by  its  long-time 
average.  If  this  approach  may  be  shown  to  be  valid,  then  a  theory  of  cumulative 
damage  which  has  been  demonstrated  to  be  valid  for  the  symmetrical  stress  pat¬ 
terns  of  the  laboratory  may  be  applied  to  the  asymmetric  stress  patterns  of  exper¬ 
ience  . 


Schjelderup  (Reference  6)  observed  further  that,  for  the  particular  stress 
records  he  examined,  there  was  considerable  consistency  in  the  normalized  prob¬ 
ability  distribution  of  alternating  stresses.  If  such  consistency  could  be  shown  to 
be  general,  tlien  it  would  be  appropriate  to  concentrate  the  effort  presently  being 
devoted  to  multi-level  fatigue  testing  into  tests  using  the  distribution  of  nature, 

The  results  of  that  type  of  testing  could  be  random  S-N  curves  which  could  be  used 
directly  for  predicting  life,  without  recovirsc  to  theories  of  cumulative  damage. 

If  Schjelderup' 8  suggestions  and  tentative  conclusions  could  be  shown  to  be 
valid,  the  number  of  parameters  which  are  required  for  describing  the  response  of 
structures  and  which  must  be  simulated  in  the  laboratory  would  be  significantly  re¬ 
duced.  These  suggestions  and  conclusions  are  therefore  examined  herein  with  the 
objective  of  providing  guidance  for  future  fatigue  analysis  and  testing. 


2.  0  RESPONSE  OF  STRUCTURE  TO  ACOUSTIC  LOADING 
2.  1  Analytic  Determination  of  Probability  Distributions 

Oscillatory  stresses  which  can  cause  fatigue  may  be  excited  by  a  variety  of 
sources.  The  sonic  fatigue  which  is  the  subject  of  this  report  is  most  frequently 
excited  by  turbulence  in  the  wake  of  a  jet  engine,  or  by  boundary  layer  turbulence. 
Both  of  these  sources  are  random  and  broad  band,  and  the  stress  response  to  those 
sources  may  therefore  be  described  most  generally  by  its  power  spectrum. 

Inspection  of  typical  spectral  densities  of  measured  or  calculated  stress  usually 
shows  the  effects  of  a  structural  filtering  process  which  transforms  the  broad  band 
excitation  into  stress  response  contained  primarily  in  a  small  number  of  rather 


Instantaneous  peak  alternating  stress  is  defined  as  1/2  the  difference  between 
a  stress  crest  and  the  succeeding  stress  trough. 

Instantaneous  mean  stress  is  defined  as  the  stress  midway  between  a  stress 
cres*^  and  the  succeeding  stress  trough. 
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narrow  frequency  bands.  The  instantaneous  stress  may  therefore  be  described 
approximately  by  Equation  (1) 


S(t)  =  Z  An(t)  Cos  (wnt +  '#>n) 


(1) 


In  Equation  (1),  Wn  is  the  natural  frequency  in  the  hfh  structural  mode.  An(t) 
is  the  instantaneous  value  of  a  random  function  with  root  mean  square  (rms)  ar^  , 
and  (f>i\  is  the  instantaneous  value  of  a  random  phase  angle.  If  the  excitation  is 
Gaussian  (normal)  and  stationary,  and  the  structure  responds  linearly,  then  it  has 
been  shown  (Reference  7)  that  An  and  <^n  vary  slowly  as  compared  to  Cos  t  , 
and  the  probability  density  of  An  is  Rayleigh  (Equation  2). 


p(An) 


(2) 


Occasionally  a  single  mode,  the  ith  >  is  excited  to  the  effective  exclusion  of 
other  natural  modes  of  the  structure.  In  that  event.  Equation  (1)  is  reduced  to 
Equation  (3) 


S  -  AiCos(cujt  +  <^|) 


(3) 


The  slow  variation  of  Aj  and  <f}\  as  compared  to  tlie  variation  of  Cos  Wit  , 
allows  Equation  (3)  to  be  construed  as  a  quasi-sinusoid,  with  negligible’’'  instanta¬ 
neous  mean  stress  and  with  peak  values  of  alternating  stress  having  tlie  same  Ray¬ 
leigh  probability  distribution  as  does  Ai  .  This  limiting  case  of  structural  re¬ 
sponse  to  acoustic  noise  is  of  particular  interest  because  most  quantitative  analyses 
of  acoustic  fatigvie  in  the  past  have  assumed  symmetrical  stress  with  Rayleigh  distri¬ 
bution  of  peaks,  and  have,  therefore,  tacitly  assumed  response  in  a  single  mode. 

More  general  and  more  realistic  response  spectral  densities,  describing  stress 
in  several  modes,  may  now  be  considered.  Two  limiting,  but  common,  cases  of 
multi-modal  response  will  be  studied  separately;  one  where  the  modal  frequencies 
are  widely  separated  and  the  other  where  the  frequencies  are  closely  spaced. 

Considering  first  the  case  where  a  number  of  well  separated  modes,  up  to  the 
K(h  ,  participate  in  the  response  to  the  same  order  of  magnitude  and  to  the  approxi¬ 
mate  exclusion  of  other  modes,  Equation  (1)  becomes  Equation  (4), 


A  steady  force  may  exist  in  addition  to  the  random  noise  excitation.  In  that 
event,  the  instantaneous  mean  stress  may  not  be  small,  but  its  variation  from 
the  average  mean  stress  will  be  small. 
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S  =  A|  Cos  (uj|t  +  <^|)  + 


^  AkCos(uk^  * 


(4) 


Because  the  An's  and  ^n's  vary  slowly,  the  time  derivatives  o£  S(t)  may  be  ap¬ 
proximated  by  Equations  (5)  and  (6), 


=  -  cj|  A|  Sin  ( (jL>|t  +  (^|)  + - ~W((AKSin(ui|<l+<^n)  (5) 

2 

H  ^  2  2 

2  ““uj|  A|Cos(Gj|ti'f^i)'i  Aj^Cos  (ai|^  t  ^  (^ ) 

dt 


Since  lias  been  assumed  (well  separated  modes)  to  be  much  greater 

than  one,  and  eacli  mode  has  been  assumed  to  be  participating  to  the  same  order  of 
magnitude,  the  rate  of  change  of  slope  (Equation  6)  is  dominated  by  the  Kjj,  term. 
Therefore,  d^S/dt^is  large  when  the  term  in  Equation  (5)  is  zero.  Consequently 

dS/dt  (Equation  5)  will  be  zero  (and  define  an  extrema  in  Equation  4)  within  a 
short  time  range  about  the  time  when  the  Ktt,  term  of  Equation  (5)  is  zero.  Thus, 
extrema  in  the  composite  stress  occur  at  time  intervals  of,  on  the  average, 

With  that  knowledge*,  tlu*  alternating  and  mean  stress  may  be  defined  by  Equations 
(7)  and  (8)  respectively; 


Sa  ® 


+  -^  2  An|Cos(WnT  +  -Cos(u;nT  +  ~  TT  +(^j|  (7) 

2  An<|Cos((i;nT +(^n)  +  Cos(uJ„T  +  TT  +  ()!>n)|  (8) 


where  T  is  any  time  when, 


oj^  T  +  =  rriTT  m  =  0, 1,  2 


(9) 


By  trigonometric  identities,  Equations  (7)  and  (8)  may  be  transformed  to 
Equations  (10)  and  (11)  respectively; 


) 

2 


ns| 


An  Sin 


Sin(iUn  T  +  (p^) 


(10) 
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S„=  i  AnCos(^Y)  Sln{w„T +<^;) 

ri8|  K 


(11) 


The  alternating  stress  given  by  Equation  (10)  consists  of  a  term,  Ak  >  with 
Rayleigh  probability  distribution,  plus  the  sum  of  (k  -  1)  terms  whose  probability 
distribution  is  determined  by  the  Rayleigh  modulation  of  a  sinusoid.  In  a  similar 
manner,  Equation  (11)  defines  the  mean  stress  as  the  sum  of  (k  -  1)  terms,  each 
term  having  distribution  determined  by  the  Rayleigh  modulation  of  a  sinusoid. 

A  proof  is  developed  in  Appendix  I  of  this  report  which  shows  that  a  Rayleigh 
modulated  sinusoid  has  normal  probability.  Because  the  sum  of  any  number  of 
terms,  each  having  normal  distribution,  has  normal  distribution,  (Reference  8) 
the  mean  stress  has  normal  probability.  The  rms  of  the  instantaneous  mean  stress, 
,  is  the  root  of  the  sum  of  the  squares  of  its  components,  as  given  by  Equa¬ 
tion  (12), 


<r 


m 


(c7-„  COS^ 


(12) 


As  in  the  previous,  single  mode  case,  the  average  vrlue  of  the  mean  stress  is 
zero  unless  a  steady-state  force  exists  in  addition  to  the  random  excitation. 

Similarly,  the  alternating  stress  is  composed  of  a  term  with  Rayleigh  distri¬ 
bution,  rms  -  ,  and  a  term  with  normal  distribution  rms  -  cTq  given  by  Equa¬ 

tion  (13), 


CTr,  = 


(o-.  Sin 


(Uk  2/ 


(13) 


The  probability  density  resulting  from  the  sum  of  a  term  with  Rayleigh  distri¬ 
bution  and  a  term  with  normal  distribution  is  developed  in  Appendix  1,  and  is  given 
by  Equation  (14), 


CQexp 


P(Sa)  = 


Zo-f 


Saotk  exp 


Z((rK^  +  (rQ^)i 


(o-K^  +  CTq^  ) 


2^  2.3/2 

2(<rx  t  CTq  ) 


+  erf 


SaO'k 


arQ^Z{(rn^  +  crQ)_ 


(14) 


Equation  (14)  is  plotted  as  Figure  2,  for  several  values  of  0'(</a'Q  .  The  compound 
distribution  is  seen  to  approach  a  Rayleigh  density  rapidly  for  the  larger  values 
would  result  from  the  assumed  large  values  of  . 

On  the  other  hand,  if  is  not  significantly  greater  than  one,  then  an  es¬ 

sential  assumption  in  the  preceding  development  is  violated  and  new  equations  must 
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FIGURE  2  COMPOUND  PROBABILITY  OF  K  +  Q 

WHEN 

K  HAS  RAYLEIGH  PROBABILITY,  rms  =  (r,^  Q  HAS  NORMAL  PROBABILITY  ,  r ms  » 


WADD  TR  61-187 


-7- 


be  derived.  The  probability  distribution  of  mean  and  alternating  stresses  for  two 
closely  spaced  modes  of  response,  the  and  kjh  ,  is  now  determined  by  a  se¬ 

parate  development.  It  will  be  convenient  in  this  development  to  use  the  functions 
(12-8),  (12+8)  ,  instead  of  Wj  ,  Wk  .  Equation  (1)  particularized  to  two  modes, 
may  be  written  in  the  new  nomenclature  as  Equation  (15). 


S  =  A)  Cos[(ii-8)t +  <P|]  +  A|,Co8[(fl  +  8)t  +  <p J 


(16) 


The  use  of  trigometric  identities  allows  Equation  (15)  to  be  rewritten  as  Equa¬ 
tion  (16), 


S 


+  A(f  +  2AjA^  Cos(2St»ipK-<^j) 


Sin(i2t  + 


with  the  phase  angle,  ,  defined: 


\lf  =  Tan 


A(,Sin(St  +  <pj)  -AjSin{St'-^i,) 
A^Cos(St♦l^J)  +  AjCos(St  -^1,) 


(16) 


(17) 


The  phase  angle  may  be  differentiated  to  show  that  jf/flis  of  the  same  order 
of  magnitude  as  8/12  ,  and  is  therefore  small  if  the  modes  are  closely  spaced. 

Consequently,  Equation  (15)  may  be  construed  to  describe  a  quasi  -  sinusoid  with 
frequency. 


a 


Q.  + 


dt  “  ■  2 


(18) 


modulated  by  the  radical  of  Equation  (16).  The  probability  distribution  of  the  radi¬ 
cal  is  the  probability  distribution  of  alternating  stress  of  the  composite  wave;  the 
mean  stress,  as  in  the  case  of  single  mode  response,  has  slight  variation  from  its 
average  value. 

The  probability  distribution  of  the  radical  is  examined  in  the  Appendix  where 
it  is  shown  to  be  Rayleigh  for  the  case  when  the  root  mean  square  of  the  jth  mode 
and  the  k}|,  mode  are  equal.  The  general  case,  presents  difficulties  in 

obtaining  a  closed  form,  or  rigorous  proof  that  the  same  distribution  holds,  but  it 
appears  likely  that  the  probability  distribution  is  Rayleigh  for  any  ratio  of  . 

The  preceeding  analytic  developments  have  determined  the  probability  distri¬ 
butions  of  mean  and  alternating  stress  for  several  idealized  cases  of  structural  re¬ 
sponse  to  random  noise.  If  these  results  are  extrapolated  to  the  general  case  of 
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structural  response,  it  may  be  concluded  that  the  probability  distribution  of  instan¬ 
taneous  mean  stress  will,  in  general,  be  normal,  and  the  distribution  of  alternating 
stresses  will  be  between  normal  and  Rayleigh,  moat  frequently  approaching  more 
closely  to  the  Rayleigh  distribution. 


2.2  Experimental  Probability  Distribution  Studies 

Data  Selection 


In  order  to  verify  the  theoretical  results  presented  above  and  to  determine  the 
effects  of  the  non-linearities  which  exist  in  real  structure,  oscillograms  of  strain 
recorded  on  typical  structure  subjected  to  high  intensity  acoustic  loading  were  ob¬ 
tained.  Oscillograms  were  obtained  from; 

1.  Douglas  Aircraft  Company,  Santa  Monica  Division 

2.  Boeing  Airplane  Company,  Wichita  Division 

3.  Convair  Division  of  General  Dynamics,  Fort  Worth 

4.  Lockheed  Aircraft  Company,  California  Division 

The  oscillograms  obtained  from  Douglas  were  recorded  during  ground  runs 
of  an  RB-66.  Boeing  provided  records  made  during  runs  of  a  B-SZ,  and  Convair 
forwarded  results  of  some  B-58  engine  runs.  The  Locklieed  records  were  taken 
from  r.omi!  specially  prepared  test  panels  made  to  simulate  typical  aircraft  structure. 
The  collected  records  represent  strain  recorded  on  Rib-Skin,  Stiffened  Skin,  Honey¬ 
comb  and  Corrugated  Web  types  of  structure,  and  therefore  include  the  major  types 
of  structure  which  luive  been  damaged  by  acoustic  noise. 

Examinations  of  the  power  spectral  density  of  strain  for  each  of  the  collected 
records  showed  that  the  records  encompassed  a  broad  range  of  modal  participation, 
from  response  in  a  single  lightly  damped  mode  to  wide  band  multi-mode  response. 
Therefore,  it  was  possible  to  select  a  total  of  24*''  records  and  cover,  with  some  dup¬ 
lication,  all  of  the  typical  structural  types  plus  a’  wide  range  of  typical  modal  response. 
Table  I  lists  the  selected  traces  and  presents  brief  descriptions  of  the  structural  type 
and  the  modal  response  represented  by  each  record. 


Data  Processing  and  Analysis 

The  selected  records  were  prepared  for  data  processing  and  analysis  by 
marking  all  of  the  significant  crests  and  troughs.  A  significant  crest  (or  trough) 
was  differentiated  from  a  minor  perturbation,  for  purposes  of  data  analysis,  by 
applictition  of  the  requirement  that  an  extremum  is  included  only  if  it  defined  an 
alternating  stress  peak  equal  or  greater  than  one-fourth  (1/4)  the  maximum  alter 
nating  stress  determined  from  a  cursory  inspection  of  an  entire  record.  Justifi¬ 
cation  of  this  truncation  procedure  is  contained  in  Appendix  II. 

The  data  contained  on  the  annotated  records,  consisting  of  approximately 
2,  000  valid  crests  and  troughs  for  each  record,  were  then  translated  into  the 


23  of  the  selected  records  represented  raw  data  and  one  represented  filtered 
data  v/hich  was  analyzed  for  a  check  on  procedures  of  data  analysis. 
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SUMMARY  OF  PROCSSS&O  STRAIN  RECORDS 
TABLE  I 


Record 

Number 

Recorded 

At 

Type  Construction 

Description  of  Modal  Response  of  Strain 

1 

Douglas 

Rib  Skin,  Aluminum 

Principally  veaponie  in  a  elngle  narrow  mode,  frequency  490 

cpio 

2 

Douglas 

Rib  Skin,  Aluminum 

Response  at  450  cps.  and  490  cps.  Root  mean  square  of 
hlghsr  frequency  mode  approximately  1,5  time,  greater  than 
root  mean  square  stress  of  lower  mode. 

3 

Douglaa 

Rib  Skin,  Aluminum 

Response  at  465  cps.  and  515  cps.  Root  mean  square  of 
higher  frequsney  mode  approximately  1,8  times  greater  than 
that  of  low  frequency  mode. 

4 

Douglas 

Rib  Skin,  Aluminum 

Response  at  485  cps,,  528  cps,  and  575  cps.  Root  mean 
square  of  lower  frequency  mode  approximately  1.3  times  that 
of  each  of  the  two  higher  frequency  modes. 

S 

Douglas 

Rib  Skin,  Aluminum 

Response  at  425  cps. ,  480  cps.  and  575  cps.  Root  mean 
square  of  center  mode  approximately  2  times  greater  than 
that  of  each  of  the  side-bands. 

6 

Douglaa 

Rib  Skin,  Aluminum 

Multi-modal  response  over  the  frequency  range  350  cps.  - 
600  cps. 

7 

Convair 

Honey  comb,  stain¬ 
less  steel 

Multi-modal  response  between  200-500  cps, ,  fairly  promi¬ 
nent  50  cps.  response. 

8 

Convalr 

Corrugatiul  Web- 
Capi  Aluminum 

Broad-band  response  between  275-260  cps.  and  fairly  promi¬ 
nent  SO  cps.  response. 

9 

Coiivalr 

CarruHatod  Web- 
Cap,  Aluiiilnuiii 

Response  in  a  moderately  narrow  band  centered  at  110  cps, 
and  nearly  wliltu  noise  below  100  cps. 

10 

Convalr 

Corriiijntc<l  Web- 
Cap,  Aluminum 

Response  In  n  brosd-band  centered  at  175  cps.  and  In  a 
narrower  band  centered  at  450  cps.  White  noise  below  100 
cps, 

11 

Convalr 

Same  as  12 

Same  as  12,  but  passed  through  a  10%  band  width  filter,  cen¬ 
ter  frequency  290  cpa. 

12 

Convalr 

Bonded-beaded 
panel.  Aluminum 

Response  primarily  in  a  narrow  band  centered  at  290  cps. 
and  a  broad-band  centered  at  50  cps. 

13 

Boeing 

Honeycomb, 

Aluminum 

Response  prlmsrily  in  s  broid-band  between  IBO  cps.  and 

400  cps.  Lower  level  multi-modal  (nearly  white)  response 
below  100  cps. 

14 

Seeing 

Honeycomb, 

Aluminum 

Multi-modal  response  from  20  cps.  to  400  cps.  Fairly 
prominent  mode  at  10  cps. 

15 

Boeing 

Honeycomb, 

Alumlnuni 

Significant  response  at  10  cps,  and  275  cps. 

16 

Boeing 

Honeycomb, 

Alumlmirtt 

Response  primarily  In  moderately  wide  band  cantered  at 

200  cps. 

17 

Boeing 

Honeycomb, 

Aluminum 

Response  In  moderately  wide  bands,  at  140  cps.  and  280 
cps.  Lower  multi-mode  response  below  100  cps. 

18 

Boeing 

Honeycomb, 

Alumlmun 

Broad-band  response  at  10  cps.  and  narrow-band  response 
at  500  cps. 

19 

Lockheed 

Rib-Skin  Test 

Panel,  Aluminum 

Narrow  modes  at  125  cps.  and  175  cps. 

20 

Lockheed 

Rlb-Skln  Test 

Panel,  Aluminum 

Narrow  bands  at  130  cps.  and  190  cps. 

21 

Lockheed 

FUed  Edge  Panel, 
Aluminum 

Moderately  broad-band  response  centered  at  120  cps. 

22 

Lockheed 

Fixed  Edge  Panel, 
Aluminum 

Similar  to  (21)  but  higher  magnitude. 

23 

Lockheed 

Fixed  Edge  Panel, 
Aliunlnum 

Narrow  band  response  at  175  cpa.  Additional  lower  level 
response  at  275  cps,  and  410  cps. 

24 

Lockheed 

Fixed  Edge  Panel, 
Aluminum 

Narrow  modes  at  150  cps.  and  270  cps. 
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punched  card  form  suitable  for  digital  computer  analyses,  A  Telemeter  Magnetics 
Inc.  "Telereader"  and  an  IBM  523  Summary  Key  Punch  were  used  for  that  opera¬ 
tion. 

A  Burroughs  220  digital  computer  was  programmed  to  receive  the  cards  con¬ 
taining  the  data  and  to  determine  for  each  record  the  average  mean,  and  the  mean 
square  of  the  crests,  troughs,  alternating  stresses,  and  mean  stresses.  Proba¬ 
bility  histograms  (using  intervals  of  1%  of  the  respective  ranges)  were  computed 
for  the  mean  and  .alternating  stresses  and  the  stress  crests  and  troughs. 


Results  and  Discussion 

The  pertinent  results  of  the  data  analysis  are  presented  in  Table  II  and 
Figures  3-15. 

Table  II,  "Some  Statistics  of  the  Strain  Records,  "  contains  the  measured 
and  calculated  values,  (in  arbitrary  units)  of  the  root  mean  square  total  stress, 

O"'  ,  alternating  stress  and  mean  stress  .  The  rms  total  stress  is 

calculated  from  the  sum  of  the  squares  of  the  stress  crests  and  troughs,  divided 
by  "{I  in  order  to  approximate  the  rms  of  the  instantaneous  total  stress,  as  is 
common  practice,  instead  of  the  rms  of  the  peak  values  of  total  stress.  The  di¬ 
rectly  computed  rms  alternating  stress  is  similarly  divided  by  to  yield  Cg  , 

When  CTo  is  normaiissed  with  respect  t.n  o  '  ,  excellenl  consistency  among 

the  records  is  observed.  Therefore,  Table  II  may  ht'  arrepted  as  providing  em¬ 
pirical  evidt'iice  that  the  rms  alternating  stress  may  be  estimated  from  the  rms 
total  stress  with  good  accuracy  and  without  the  necessity  for  performing  the  formal 
procedure  of  resolving  total  stress  into  alternating  and  mean  components. 

The  normaliisation  does  not  reveal  a  similar  quantitative  consistency  in  the 
mean  stress  component.  It  does,  however,  confirm  that  for  the  records  examined, 
the  variance  of  the  mean  stress  is  substantially  smaller  than  the  variance  of  the 
total  stress. 

Table  11  also  contains  the  range  of  stress  crests  and  stress  troughs  and  the 
ratio  of  the  ranges.  These  values  are  included  because  the  ratio  of  the  ranges  is 
a  convenient  measure  of  the  degree  of  one  type  of  non-linearity  which  occasionally 
exists;  ratios  substantially  different  from  one  imply  considerable  difference  in 
stiffness  for  positive  stress  and  for  negative  stress.  It  is  observed  that  the  non- 
linearities  have  no  apparent  effect  upon  the  normalized  alternating  stress. 

Figures  3  to  9  present  the  measured  probability  integrals  of  (normalized) 
alternating  stresses  plotted  on  truncated  Rayleigh  probability  paper,  so  that  a 
measured  distribution  which  conforms  with  the  theoretical  probability  distribution 
will  plot  as  a  straight  line  (see  Appendix  II).  The  conformity  of  the  test  data  to 
straight  lines  is,  in  general,  very  good.  The  greatest  deviations  from  straight 
lines  occur  near  the  truncation  level  of  strain  and  the  high  value  (low  probability) 
levels  of  strain.  Deviations  from  linearity  near  the  truncation  can  easily  be  attri¬ 
buted  to  human  errors  in  data  processing,  which  result  from  the  euinotator  of  the 
record  attempting  to  be  conservative  and  including  some  levels  which  were  slightly 
below  the  nominal  truncation  level.  Data  points  representing  high  stress  levels 
have  an  apparent  tendency  to  fall  below  the  line.  If  that  tendency  is  real,  evidence 
of  peak  clipping  is  presented;  on  the  other  hand,  the  high  stress  points  are  derived 


WADD  TR  61-187 


-11- 


1 

1 


I- 

I 


TABLE  11 

Some  Statistics  of  the  Strain  Records 


Record 

Number 

or' 

RMS 

Total 

Stress 

CTa 

RMS 

Alt. 

Stress 

m 

O-m 

RMS 

Mean 

Stress 

a 

Rc 

Range 

of 

Crests 

Rt 

Range 

of 

Troughs 

m 

1 

63.  0 

61.5 

.98 

22.5 

.  36 

241 

208 

1.  15 

2 

65.5 

62.5 

.96 

32.2 

.49 

286 

280 

1. 02 

3 

58.0 

56.  0 

.97 

21.2 

.  37 

281 

232 

1.21 

4 

186.0 

.98 

54.  0 

.29 

619 

669 

.92 

5 

.97 

48.4 

.30 

585 

571 

1. 02 

6 

87.5 

82.5 

.95 

41.8 

.48 

387 

390 

.99 

7 

176.0 

174.0 

.  99 

53. 0 

.  30 

632 

617 

1.03 

8 

144.0 

138.0 

jm 

59.  0 

.41 

439 

813 

.  54 

9 

168.0 

162.0 

Hi 

65.0 

.39 

484 

836 

.58 

10 

283.0 

272.0 

.96 

114,0 

.40 

1044 

1103 

mm 

12 

324.0 

319.0 

.98 

89.0 

.27 

855 

934 

.  92 

13 

469.0 

rnim 

.98 

144.0 

.31 

1387 

1486 

.90 

14 

113.0 

Warn 

.95 

65.0 

.57 

500 

598 

.  84 

15 

314.0 

304.  0 

.97 

108.  0 

.  34 

1102 

879 

1.25 

16 

330.0 

326.0 

.99 

82.0 

.25 

233 

1112 

1.20 

17 

506.0 

490.  0 

.97 

179.0 

.  37 

1965 

1245 

1. 57 

18 

198.  0 

194.  0 

.98 

61.5 

.31 

673 

611 

1.  10 

19 

268.  0 

264.0 

.98 

64.  5 

.24 

748 

939 

00 

o 

20 

213.0 

210.0 

.99 

53.0 

.  25 

564 

646 

.86 

21 

159.  0 

154.  0 

.97 

55.5 

.  35 

486 

662 

.86 

22 

229.0 

218.0 

.95 

97.  0 

796 

955 

.83 

23 

182.  0 

178.0 

.98 

52.5 

.29 

590 

661 

1. 06 

24 

183.0 

179.0 

.98 

54.  5 

.  30 

645 

556 

1.  16 

Average 

Standard  Deviation 

mm 

Hlliii 
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on  tihe  basia  of  a  small  number  (2-lOy  of  measurements  and  the  confidence  which 
may  be  attached  to  their  location  is  doubtful. 

It  is  especially  encouraging  to  observe  the  probability  distribution  for 
records  8  emd  9  in  Figures  4d  and  5a.  Those  records  represent  stress  response 
which  is  quite  non-linear,  as  shown  by  ratio  of  range  of  crests  to  reinge  of  troughs 
in  Table  II.  However,  the  probability  integrals  for  those  records  conform  to 
straight  lines  at  least  as  well  as  does  the  average,  demonstrating  that  the  effects 
of  that  particular  type  of  non-linearity  on  the  probability  distribution  of  alternating 
stresses  is  small. 

Figures  7d  and  8a  for  records  21  and  22  show  considerable  deviation  from 
the  straight  line.  Those  records  were  obtained  from  Lockheed  and  pertain  to 
measurements  made  on  speciall-y  constructed  panels  rather  than  on  actual  aircraft 
structure.  No  further  explanation  for  the  deviation  of  the  measured  results  from 
the  theoretical  results  is  at  hand. 

Figure  9  is  the  superposition  of  Figures  3-8  (with  the  exception  of  Figure  7d 
and  8a  where  some  unknown  factors  were  apparently  operating).  This  figure  sum¬ 
marizes  the  great  consistency  of  the  probability  distributions  and  the  normalized 
standard  deviation  of  alternatiihg  stress,  as  presented  in  the  preceding  figures  and 
in  Table  II.  The  lack  of  scatter  contained  in  Figure  9  should  be  considered  in  light 
of  the  variety  of  structural  configurations  and  modal  response  represented  by  the 
records  it  summarizes.  The  consistency  in  spite  of  those  differences  is  presented 
fis  sufficient  evidence  that  the  normalized  level  and  probability  distribution  of  alter¬ 
nating  streaao.s  are  consistent  factors  in  aircraft  subjected  to  high  intensity  acous¬ 
tic  noise. 

Figures  10-16  present  the  measured  probability  integral  of  mean  stress  on 
normal  distribution  paper.  Again,  the  correspondence  to  straight  lines  is  good. 
Figures  lid  and  12a  show  significant  deviation  from  straight  line.s.  Those  figures 
pertain,  as  did  Figures  4d  and  5a  to  the  distinctly  non-linear  records,  8  and  9.  It 
may  be  useful  to  observe  that  this  type  of  non-linearity  has  more  apparent  effect 
upon  the  mean  stretci  than  on  the  alternating  stress. 

No  composite  plot  of  mean  stress  probability  distribution  is  presented.  Such 
a  plot,  if  made,  would  reflect  the  lack  of  uniformity  of  the  normalized  rms  mean 
stress  tabulated  in  Table  II  und  would  show  considerable  scatter. 


3.0  APPLICATIONS 

It  appears  to  be  established  by  theory  and  experiment  that  it  is  possible  to 
describe  the  probability  distribution  of  alternating  stresses  of  structures  subjected 
to  high  intensity  acoustic  loading  by  a  simple,  consistent,  and  conveniently  repro¬ 
duced  function.  This  could  have  great  influence  in  optimizing  the  techniques  of 
future  fatigue  testing. 


3.0.1  Qualification 


Before  examining  the  ramifications  of  that  conclusion,  it  is  appropriate  to 
emphasize  that  a  procedure  which  resolves  the  total  stress  into  its  mean  stress 
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component  and  alternating  stress  component  is  Justified  only  if  it  can  be  established 
that  but  one  of  the  components  is  significant  in  fatigue.  This  has  not  yet  been  estab¬ 
lished.  Consequently,  the  usefxilness  of  the  results  presented  in  Section  U  above 
and  the  vcilidity  of  the  applications  which  will  be  presented  below  remain  subject  to 
a  test  program  which  should  be  initiated  as  soon  as  possible.  Appendix  III  of  this 
report  outlines  test  procedures  which  will  be  sufficient  for  examining  the  assump¬ 
tion  that  the  instantaneous  value  of  the  varying  mean  stress  has  only  second  order 
effects  upon  the  life  of  a  structure  subjected  to  random  stresses.  In  the  remainder 
of  this  section  it  is  assumed  that  the  required  test  program  has  been  performed 
with  favorable  results,  and  the  presence  as  well  as  the  effects  of  the  varying  mean 
stress  may  therefore  be  neglected. 


3.  1  Life  Prediction 

Theories  of  cumulative  damage  which  are  more  sophisticated  than  the  Palmgren 
Miner  Hypothesis  (Reference  1),  have  been  proposed  by  several  independent  re¬ 
searchers  (Re  ferences  2,  3,  4).  Although  those  theories  make  use  of  apparently 
markedly  different  assumptions  concerning  the  mechanism  of  interaction  between 
stresses  at  different  levels,  they  (and  others)  yield  remarkably  similar  results. 
Essentially,  modern  theories  of  cumulative  damage  agree  that  the  existence  of  a 
range  of  stresses  has  the  qualitative  effect  of  causing  "rotation"  of  the  basic  S-N 
curve  of  a  material.  If  the  amount  of  "rotation"  were  known  quantitively,  then 
the  elementary  theory  could  be  applied,  using  the  rotated  S-N  curve,  for  suitably 
accurate  random  fatigue  life  prediclion.s.  Much  of  the  effort  prebently  being  expen¬ 
ded  on  the  development  of  cumulative  damage  theories  is  focussed  upon  deriving 
techniques  for  predicting  the  degree  of  rotation  caused  by  stress  interaction.  The 
experimental  methods,  in  all  cases,  make  use  of  tests  wherein  two  or  more  levels 
of  stress  are  applied  in  random  sequence  to  test  specimens,  The  probability  distri¬ 
butions  of  the  applied  stresses  arc  often  arbitrary,  or  are  intended  to  simulate  the 
probability  distribution  of  gust  loading. 

It  may  be  postulated  that  the  empirical  expression  for  rotation  of  the  basic 
S-N  curve  which  is  derived  from  results  of  tests  using  a  particular  probability 
distribution  of  stresses,  will  provide  the  best  possible  accuracy  for  cases  where 
the  stress  probability  distribution  of  usage  is  identical  to  the  probability  distribution 
of  the  lest.  It  follows  from  that  postulate,  and  because  the  probability  distribution 
of  (alternating)  stresses  in  life  can  be  described  with  engineering  accuracy  as  a 
Rayleigh  distribution,  that  testing  performed  in  the  future  with  the  intent  of  deri¬ 
ving  a  theory  of  cumulative  damage  for  use  in  acoustic  fatigue  should  incorporate 
a  Rayleigh  distribution  of  stresses. 

The  results  (life  versus  rms  stress)  of  tests  performed  with  the  Rayleigh 
probability  distribution  of  stresses  may,  in  addition  to  being  incorporated  into  a 
theory  of  cumulative  damage,  be  used  directly.  Specifically,  if  a  random  S-N 
curve  for  a  particular  probability  distribution  is  available  and  there  is  reason  to 
believe  that  a  structure  will  be  subjected  to  stresses  with  a  similar  probability 
distribution,  then  there  is  no  need  for  a  theory  of  cumulative  damage;  it  is  only 
necessary  to  determine  the  rms  alternating  stress  and  the  effective  frequency  and 
the  life  may  be  read  directly  from  the  random  S-N  curve.  That  procedure  is 
obviously  simpler  for  designers  to  use  than  a  theory  of  cumulative  damage  can 
reasonably  be . 
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3.2  Production  of  a  Rayletgh  Distribution 
3,2.1  Production  of  a  Continuous  Distribution 

A  benefit  accruing  from  the  fact  that  the  probability  distribution  of  (alter¬ 
nating)  stresses  has  been  found  to  be  Rayleigh  is  the  ease,  in  many  cases,  of  re¬ 
producing  that  distribution.  It  may  be  produced  in  the  laboratory  by  taking  advan¬ 
tage  of  the  same  natural  process  which  produces  it  in  actual  structure:  that  is, 
mechanical  filtering  of  a  broad  band  excitation.  The  natural  frequencies  of  a  test 
specimen  may  be  determined,  and  the  specimen  may  then  be  subjected  to  random 
vibration  with  a  frequency  content  centered  approximately  at  the  specimen  natural 
frequency  which  it  is  desired  to  excite,  and  with  a  band  width  substantially  greater 
than  the  anticipated  specimen  bandwidth.  The  response  stress  of  the  specimen  will 
then  approach  a  perfect  Rayleigh  distribution. 


3.2,2  Simulation  by  a  Discrete  Distribution 

In  many  laboratories,  test  equipment  limitations  will  preclude  the  use  of 
resonant  response  to  random  excitation  for  production  of  a  Rayleigh  distribution. 

In  that  event,  it  is  desirable  to  define  a  test  performed  at  a  finite  number  of  dis¬ 
crete  stress  levels  wliich  will  produce  the  same  effects  as  a  Rayleigh  distribution. 
That  problem  would  become  trivial  if  an  adequate  theory  of  cumulative  damage  were 
in  existence,  since  a  straightforward  application  of  a  law  of  cumulative  damage 
could  be  the  calculation  of  an  cfiuivalent  stress,  Se  such  that  application  of  N 
cycles  of  random  stress,  rrns  ~  cy  .  or  N  cycles  of  constant  stress,  Se  ,  would 
cause  failure.  In  reviewing  existing  theories  of  cumulative  damage,  it  is  natural, 
because  of  the  ease  of  application,  to  investigate  first  the  results  of  using  the  ele¬ 
mentary  theory.  When  Miner's  Law  (Reference  1)  is  used,  the  accuracy  of  predic¬ 
tion  of  Se  ranges  from  fairly  good  to  poor;  the  goodness  of  fit  is  found  to  be,  as 
would  be  expected,  best  when  the  range  of  simulated  random  stresses  is  low  and 
the  S-N  curve  is  well  known  and  poor  v/lion  the  range  of  random  stresses  is  broad 
and/or  stress  concentrations  in  the  structure  prevent  an  adequate  knowledge  of  the 
S-N  curve, 

A  detailed  review  of  the  more  elaborate  theories  of  cumulative  damage  in 
fatigue  does  not  contain  promise  of  a  better  method  for  determining  a  single  level 
So  from  a  loading  probability  and  an  S-N  curve.  Instead,  study  of  the  reasoning 
behind  the  more  elaborate  cumulative  damage  theories  has  suggested  that  discrete 
amplitude  testing  which  will  enable  the  life  of  a  random  loaded  structure  to  be 
simulated  must  involve  more  than  one  level  of  stress,  with  the  sequence  of  levels 
mixed  in  a  random  manner.  The  highest  level  of  discrete  stress  must  be  near  the 
highest  level  of  random  stress  which  occurs  significantly  often  and  the  lowest  level 
of  discrete  stress  must  be  approxinjately  as  low  as  the  lowest  level  of  random  or 
approximately  33%  of  the  highest  stress,  whichever  is  higher.  More  than  two 
levels  of  discrete  stress  will  be  required  if  the  S-N  curve  is  not  well  known,  eitlier 
because  of  stress  concentrations  or  because  of  the  effects  of  stress  interaction. 

If  a  very  large  number  of  discrete  levels  of  stress  are  used,  the  test  would 
approach  exact  simulation  and  no  knowledge  of  the  effects  of  interaction  or  stress 
raisers  would  be  required.  However,  the  difficulties  of  performing  multi-level 
tests  would  soon  approach  and  surpass  difficulties  of  performing  true  random  tests 
and  the  desired  economic  advantages  of  discrete  testing  would  be  lost.  The  fol¬ 
lowing  is  proposed  as  a  compromise; 
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"1.  It  Is  assumed  that  the  effects  of  interaction  between  stresses  jit  level 
Smax  stresses  at  level  Smtn  are  an  order  of  magnitude  greater  than  the" 
effects  of  interaction  between  stresses  with  a  range  of  E5%  to  33%  as  great. 

That  assumption  is  in  accordance  with  phenomenologicail  theories  of  fatigue  as 
proposed  by  Freudenthal  (Reference  2),  Fuller  (Reference  3)  and  others  and  is 
in  accordance  with  considerable  experimental  data. 

2.  The  assumption  in  paragraph  1  allows  the  statement  that  an  adequate  overall 

simulation  may  be  achieved  provided  the  interaction  between  extreme  stress  levels 
is  simulated;  interaction  within  fractions  of  the  extreme  cause  second  order  effects 
which  may  be  ignored.  That  assumption  further  allows  the  statement  that  a  form 
of  Miner's  Law  is  usefully  accurate  for  predicting  the  stress,  S',  which  is  equiva¬ 
lent  to  stresses  within  the  narrow  range  where  interaction  is  not  significant.  In 
the  form  used,  1 

Sj  1'b 

/  S%lS)dS 

_Si _ 

p(S)dS 

Si 

where  p  (S)  is  the  probability  density  of  stresses  and  b  is  the  inverse  slope  of  the 
S-N  diagram  in  log-log  representation.  (If  S.  -  Sj  is  reasonably  narrow,  a  straight 
line  representation  will  be  completely  adequa'te.) 

3.  If  H  range,  Sj  -  Si  may  be  found  so  that  S'i-j  is  essentially  independent  of  b 
over  a  realistically  wide  range  of  b,  then  S'  may  be  specified  without  quai:titalive 
knowledge  of  b.  Material  variability,  stress  concentrations,  and  stress  inter¬ 
actions,  etc.  which  affect  b  thus  become  second  order  effects  in  evaluating  the  ac¬ 
curacy  of  the  life  which  is  determined  by  the  testing  procedure. 

To  that  end,  the  values  of  S'  which  are  intended  to  simulate  random  stresses 
between  Sj  and  Sj,  for  the  Rayleigh  probability  density  and  several  realistic  values 
of  b,  have  been  computed.  The  results  follow; 


TABLE  III 

Single  Stress  Level  Equivalent  to  a  Segment  of  a  Rayleigh  Density  of  Stress 


Si 

Sj 

S'{b=4) 

S’(b=6) 

S'(b=8) 

0 

cr 

.  735 

.  773 

.  80 

CT 

2  cr 

1.51 

1.  56 

1.60 

2cr 

30- 

2.  39 

2.42 

2.45 

3cr 

oo 

3.34 

3.38 

3.41 

(need  not  be  included  in 
test  procedure) 
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The  differences  in  stress  levels  within  tlie  rows  of  Table  in  ere  seen  to  be 
small,  indicating  that  an  accurate  knowledge  of  the  shape  of  the  S-N  curve  is  of 
second  order  importance  for  a  realistically  wide  range  of  b.  Consequently,  it  is 
possible  to^rescribe  a  three»levol  test  which  is  nearly  equivalent  to  a  random 
test.  Typically,  if  it  is  desired  to  simulate  a  Rayleigh  distribution,  rms  =  c  , 
it  could  be  done  by  means  of  blocks  of  stress,  each  block  consisting  of  1,110 
repetitions  of  stress  equal  to  3.4(r  ,  12,440  cycles  of  stress  equal  to  2. 4  O'  ,  and 
47,  100  cycles  of  1.6  O’  .  Each  block  would  be  equivalent  to  100,  000  cycles  of 
random  stress  (39,  350  cycles  of  stress  at  levels  less  than  <r  may  be  ignored)  and 
the  sequence  of  high,  low  and  medium  stresses  would  be  varied  within  each  succeed¬ 
ing  block  in  a  random  manner.  N  blocks  would  be  equivalent  in  life  to  100,  OOON 
cycles  of  random  stress. 


3. 3  Proof  Testing 

It  is  reasonable  to  assvime  that,  even  after  a  reliable  theory  of  cumulative 
damage  in  fatigue  is  developed  it  will  continue  to  be  appropriate  to  perform  fatigue 
tests  on  subassemblies,  major  assemblies,  and  complete  structures,  in  order  to 
evaluate  specific  designs  and  a.ssembly  techniques.  Simulation  of  the  probability 
distribution  of  stresses  is  more  important,  if  anything,  for  this  type  of  testing  than 
it  is  for  the  laboratory  type  experimentation  discussed  in  Section  3.2  above.  Simu¬ 
lation  is  both  automatic  and  perfect  when  the  proof  testing  is  accomplished  by  sim¬ 
ply  tying  down  the  tost  airplane  and  letting  the  engine  run;  however,  when  economic 
considerations  dictate  testing  of  local  areas  using  siren  excitation,  then  the  consid- 
ei'atiuiis  of  paragraph  3.2.2  apply  and  a  simulated  random  test  should  be  performed. 


3.4  Accelerated  and  Aggravated  Testing 

Fatigue  proof  testing  of  aircraft  structure  is  accelerated  testing  whenever 
it  is  assumed  that  the  stress  history  during  flight  is  ao  trivial  compared  to  the 
stress  history  during  the  take  off  run,  that  simulating  the  take  off  stress  and  the 
total  time  at  the  stress  level  is  equivalent  to  simulating  the  entire  life  of  the  air¬ 
plane.  Granting  that  assumption,  proof  testing  at  maximum  ground  levels,  for 
perhaps  one  hundred  hours,  i.s  considered  to  qualify  an  airplane  design  for  its  en¬ 
tire  life  of  up  to  10,  000  flight  hours. 

Although  an  important  amount  of  acceleration  of  airplane  life  is  contained  in 
testing  at  normal  ground  stress  levels,  there  is  a  natural  desire  to  shorten  test 
times  still  further.  This  may  be  accomplished  with  reasonable  accuracy  provided 
the  random  S-N  curve  of  the  material  or  joint  being  tested  is  known  and  it  is  pos¬ 
sible  to  apply  an  exaggerated  stress  level.  This  technique  is  discussed  at  some 
length  by  Schjelderup  (Reference  9)  and  will  not  be  considered  further  here,  except 
to  remark  that  both  normal  test  techniques  which  assume  negligible  effects  upon 
life  from  the  stress  levels  of  cruise  auid  the  aggravated  tests  described  by 
Schjelderup  (Reference  9)  assume  better  knowledge  of  random  S-N  curves  than  is 
usually  at  hand.  Substantial  improvement  in  confidence  that  may  be  attached  to 
accelerated  testing  would  result  from  a  greatly  enlarged  library  of  random  S-N 
curve  s . 
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4.  0  CONCLUSIONS  AND  RECOMMENDATIONS 


Aqv  stress  history^  no  matter  how  complex,  may  be  resolved  into  sun  alter¬ 
nating  stress  superimposed  upon  a  varying  measi  stress.  The  resolving  procedure 
is  essentially  a  high-pass  filtering  process  for  the  alternating  stress  and  a  low- 
pass  process  for  the  mean  stress. 

The  filtering  tends  to  make  the  probability  distribution  of  the  alternating 
stress  component  conform  approximately  to  the  probability  distribution  charac¬ 
teristic  of  the  highest  frequency  mode  which  is  present  to  a  significant  degree  in 
the  total  stress.  Therefore,  the  probability  distribution  of  peaks  of  the  alternating 
stress  component  approaches  a  Rayleigh  distribution  when  the  stress  is  excited  by 
acoustic  noise. 

Similarly,  the  probability  distribution  of  instantaneous  values  of  the  stress 
in  lower  frequency  mode  s  tends  to  dominate  the  distribution  of  instantaneous  tnean 
stresses.  Consequently,  a  normal  distribution  of  mean  stresses  tends  to  occur. 

Examination  of  oscillograms  containing  strain  measured  upon  a  wide  variety 
of  typical  aircraft  structures  exposed  to  high  intensity  acoustic  noise  has  confirmed 
that  the  expected  probability  distribution  of  mean  and  alternating  stresses  are  ob¬ 
served,  within  limits  of  engineering  accuracy. 

Resolving  total  stress  into  mean  and  alternating  coinjjonents  will  be  extremely 
useful  provided  it  may  bo  demonstrated  by  test  that  the  instantaneous  value  of  mean 
stress  has  only  slight  effect  vipon  the  fatigue  damage  caused  by  the  alternating  stress 
component.  A  test  program  to  determine  the  effect  of  variable  mean  stress  with 
variance  of  1%  to  30%  the  variance  of  the  alternating  stress  should  therefore  be 
initiated  in  the  immediate  future.  An  outline  of  a  sufficient  test  program  for  deter¬ 
mining  those  effects  is  included  in  tliis  report  as  Appendix  HI. 

If  the  results  of  the  recommended  test  program  are  favorable,  then  future 
fatigue  testing  of  laboratory  specimens  and  real  structures  should  reflect  the 
probability  distribution  of  which  has  been  found  to  be  nearly  universal  for  struc¬ 
tures  excited  by  high  intensity  acoustic  noise.  The  results  of  such  fatigue  testing 
will  be  random  S-N  curves  which  may  be  used  directly  in  the  design  and  testing  of 
structures  and  which  may  be  used  as  data  for  establishinga  general  theory  of  cumu¬ 
lative  damage  in  fatigue. 
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APPENDIX  1 

COMPOUND  PROBABILITIES 


In  paragraph  2,  preceding,  the  form  of  the  mean  stress  and  alternating  stress 
was  found  to  be  the  product  and/or  sum  of  two  or  more  random  and  mutually  inde¬ 
pendent  functions.  The  probability  functions  resulting  from  those  compound  results 
are  derived  herein. 


1.  Consider  the  function,  S  =  A  Cos  fl  where  A  has  Rayleigh  probability  with 
variance  of  ffg  ,  and  all  values  of  0  are  equally  probable. 

Since  all  values  of  6  are  equally  likely,  the  probability  that  Cos  9>^  is, 


cosfl 


It/  ' 

r 

-  Cos  Y 

t  Q 
-|< 

=  0 

.Since  A  lias  a  Rayleigh  density,  the  probability  that  A  is  bounded  by  T,  T  +  dT  is 


Pj(r)  =  0 


T  <  0 


=  ^  exp 


r  >0 


If  A  is  bounded  by  T,  T  +  dT  ,  and  Cos  6  ,  then  A  Cos  9  >  Q.  The  proba¬ 

bility,  then,  that  A  Cos  0  >  Q,  is  the  product  of  the  above  mutually  independent 
probabilities,  summed  over  all  values  of  T  , 


P(Q)  =  4- 


Cos 


Q 

f  e>tP 


-dr 


The  above  expression  for  the  Probability  Integral  of  A  Cos  8  ,  is  not  integrated 
easily.  However,  if  its  derivative  is  taken, 
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^  p(Q)  s  X 

dQ  TT 


Texp 

.jlL 

L^J 

2  nf 

0*0  V 

r  -Q 

IQI 


dT 


the  integration  nnay  then,  be  performed  and  multiplied  by  (-1)  to  yield  the  proba¬ 
bility  density  of  A  Cos  6  , 


which  is  recognized  as  the  normal  probability  density. 


Z.  Consider  the  function,  A  +  B,  where  A  has  Rayleigh  probability,  variance  of 
CTg^  ,  and  B  has  normal  probability,  varijmce  of  cr|j^  .  Since  A  has  Rayleigh 
distribution  the  probability  that  A  >  (Q  -  T  )  i*;  , 


P^(Q-T)  =  exp 
=  I 


(Q'T) 


2or 


0  J 


T  <  0 
T  >  0 


B  has  normal  density.  Therefore,  the  probability  that  B  is  bounded  by  TiT  +  dT  is, 


Pb^'t)  “ 


I 

o-b 


exp 


dT 


If  A  >  Q  -  T  ,  and  B  is  bounded  by  T,  r  +  dT,  then  A  +  B  >  Q,  Since  values  of  A  and 
B  are  assumed  mutually  independent,  the  joint  probability  is  the  product  of  the 
probabilities.  Summed  over  allT  ,  the  re.sultant  probability  is, 


Pa,b(Q)  ' 


The  integration  is  straightforward  and  the  result  is, 


Pa.8(Q)= 


9 

(Toexp 

Q 

2(or/  +  or/) 

2(<ra  +  CTb  ) 


1  +  erf 


l-erf 
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or,  taking  the  negative  derivative  to  yield  the  probability  density, 


Qo-jexp 


2^  2.3/2 

2(0-0  +  O'b  ) 


1  +erf 


Qo'fl 

<Tb®*P 

'  q"  ' 

a-bV2((r®  +  a^) 

3. 


Consider  the  function  A 

.2 


+ 


2AB  Cos  d 


bution,  variance  CTo*  ,  B  has  Rayleigh  distribution 
of  d  are  equally  probable. 


where  A  has  Rayleigh  distri- 
,  variance  Cb*  .  All  values 


2  2 

The  probability  integral  of  (A  +  B  )  may  be  found  readily,  as  may  the  proba¬ 
bility  density  of  2AB  Cos  6  .  However,  reflection  shows  that  it  is  not  correct  to 
perform  the  convolution  of  P(A2  +  B^)  with  p(2AB  Cos  8  )  to  determine  P(A2  + 

+2AB  Cos  6  )  because  the  instantaneous  value  of  (AB)  is  not  independent  of  the  in¬ 
stantaneous  value  of  (A2  +  b2).  That  is,  2AB  is  equal  to  or  less  than  a2  +  b2. 


It  is  necessary  to  make  use  of  trigonometric  identities  to  write, 


A^+  2AB  Cose  =  [a^+  B^_  j  t  Sin(2Tan”'  Cos0_ 


2  2 

In  that  form,  the  instantaneous  value  of  (A  (•  B  )  does  not,  in  any  way,  govern 
or  limit  the  instantaneous  value  of  A/B,  and  a  convolution  integral  may  therefore  be 
performed  on  the  (now  independent)  bracketed  terms.  That  is,  if  A^+B^  2  Q/T  , 
and  1  i  Sin  2  Tan"  _  A/B  Cos  8  is  bounded  by  T  ,  TtdT  ,  then, 

[a2  +  [l  +  SinETan-l  A/B  Cos»]  >Q. 

2  2  2  2 
The  probability  that  A  -t-  3  >  Q/T  will  first  be  determined;  A  +  B  will  be 

greater  than  Q/T  ,  provided 


A^  >  ^  -  V 


B  bounded  by  V,  V  +  dV, 


Since  A  and  B  each  have  Rayleigh  distribution,  the  above  probabilities  may 
be  written  immediately: 


WADD  TR  61-187 


34- 


P,a(^r-V)  =  .«pr- 


(Q/t-V) 


2<rfl 


»  I 
«»  0 


I 


20-fc 


exp 


20-l 


dV 


V  <  4 

V  >  4 

V  <  0 

V  >  0 


Therefore, 


1 

Q/t 

.  0  1 

eo 

dV  sxp 

q/t 

•  -1 
V 

#IV/ 

20-/ 

1  «*P 

"ii 

2Tcr/_ 

uV 

which  may  be  integrated  directly  to  yield 


CTb  exp 


2T0rL 


-OTg  exp 


2Tcr 


0  J 


2  2 
-0-0 


-1 


The  probability  that  1  +  SinZTan  A/B  Cos  6  is  bounded  by  T,T  +  dT  is  now 
found.  This  is  done  by  first  finding  the  probability  integral  of  (A/B),  determining 
from  it  the  probability  integral  of  Sin2Tan~ ^A /B,  and  then  performing  a  convolu¬ 
tion  integral  of  that  probability  with  tlie  probability  density  of  Cos  6  . 

Proceeding  as  in  the  previous  cases,  it  is  seen  that  A/B  >  Z  provided 
A  >  Z  ^  ,  auid  B  bounded  by  ^  -d<^,  ^  .  A  amd  B  have  Rayleigh  probabilities; 
therefore 


P^(Z^)  =  6xp 
PbU) 

Pa(Z)' 


2<7-® 


u 

*  -1 

2  exp 
O’b 

d^ 

r* 

Ok 

A  A 

1  i> 

I  AV  IS 

zV 

1 

>0 

>0 
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which  may  be  integrated  immediately  to  yield, 


PsuZTon”'  -ICoi^CH)  “  g*  (l" 


R(o-q^  +  tr^) 


i 


R*  {a-^  +  cr^)+  cr^ 


Now,  if  Sin2Tan"^A/B  >R,  then  1  +  SinZTan"  /BCoa  fl  :»  1  +  R. 

Define  1  +  R  aT  ,  and  it  is  seen  that  the  probability  that  1  +  Sin2Tan"^A/B 
Cos  e  >  R  +  1  is, 


p(t)  =i(i- 


(r-DCo-o^  +  a^) 


2x2  ,2  2x2,  2  ^  X 

o’b  )  (o-Q  “  o’b  )  "  2t) 


Tine  negative  differential  of  the  above  function  may  be  taken  to  yield  the  prob¬ 
ability  that  (1  +  SinZTan'^A /BCos  0  )  is  between  T,  T  +  dT  . 


P(t)  = 


-  2  2  .  2  2^ 
2o-0  (Tb  (CTfl  +  CTb  )  dt 


/,  x*  2  2  2  2  2 

(I-t)  (o-o  -o-b  )  4cr„  o-b 


3/2 


The  indicated  quadrature  may  be  performed  readily  in  closed  form,  for  the 
special  case  ctq  ■  <rb  ,  and  yields  a  Rayleigh  probability  integral, P(Q)=e;^?|jQ/2((rQ  t  o-b^)] 
Efforts  to  perform  the  integration  in  closed  form  for  the  general  case,  ffa  Cb  , 
have  not  been  successful.  However,  numerical  integration  yielded  results  for  the 
points  calculated  which  were,  again,  in  conformance  with  a  Rayleigh  distribution. 

It  is  extremely  likely  that  the  solution  to  the  general  case,  if  it  could  be  found,  has 
a  Rayleigh  distribution. 


WADD  TR  61-187 


-37- 


APPENDIX  II 

TRUNCATION  AND  RECTIFICATION  OF  DATA 


Data  Processing 

When  a  strain  record  contains  the  results  of  response  in  several  modes, 
there  may  be  a  large  number  of  points  of  zero  slope  separated  by  small  changes 
of  strain.  At  present  no  universally  accepted  criterion  exists  for  deciding  which 
of  those  extrema  are  significant  in  fatigue  analysis. 

It  is  certainly  sufficient  to  make  no  distinction  between  minor  perturbations 
and  significant  fluctuations  and  to  include  all  points  of  zero  slope  in  a  data  analysis 
of  stress  peaks.  The  question  of  whether  it  is  necessary,  as  well  as  sufficient, 
arises  when  a  high  frequency  low  amplitude  mode  is  excited  along  with  lower  fre¬ 
quency  high  amplitude  modes.  That  condition  of  modal  response,  together  witli 
the  assumptions  that  stress  peaks  are  the  damaging  mechanism  in  fatigue  and  all 
points  of  zero  slope  in  a  stress  record  are  significant  pe2dcB,  led  Weibull  (Refer¬ 
ence  5)  to  conclude  that  every  stress  record  is  an  entity  which  c  annot  be  compared 
to  anything  except  a  nearly  identical  stress  record. 

A  dliferent  difficulty  arises  when  considering  alternating  stress,  instead 
of  stress  peaks,  as  the  damaging  mechanism.  In  that  case,  inclusion  of  all  extre¬ 
ma  can  transform  a  normal  stress-time  history  into  a  superposition  of  a  small 
alternating  stress  component  upon  a  large  mean  stress  component.  That  condition 
is  the  reverse  of  what  was  considered  when  the  Alternating  Stress  Hypoth(?sis 
(Reference  6)  was  proposed  and  it  may  be  concluded  that  either  we  may  not  include 
all  minor  extrema  or  the  alternating  stress  hypothesis  must  be  abandoned. 

A  common  technique  for  discarding  minor  perturbations  involves  assuming 
that  fluctuations  which  do  not  cause  zero  crossings  are  trivial,  and  a  valid  crest 
is  the  largest  excursion  from  zero  between  two  adjacent  zero  crossings.  This  as¬ 
sumption,  however,  is  frequently  self-defeating.  For  example,  consider  the  ex¬ 
treme  case  when  a  significant  steady-state  stress  exists.  In  that  event,  many 
(or  all)of  the  stress  fluctuations  will  not  cross  true  zero  and  will  be  ignored.  That 
case  may  perhaps  be  handled  by  (mentally)  filtering  out  the  steady  component 
and  sketching  in  aline  which  apparently  represents  the  average  stress.  Using 
that  line  as  a  first  approximation  to  zero,  instantaneous  mean  stress  (point  mid¬ 
way  between  a  "line  defined"  crest  and  a  "line  defined"  tre  .gh)  may  be  tabulated 
and  a  second  approximation  to  zero  may  be  determined  by. 


N 

Z  Sm; 

M 

N 


However,  the  process  may  converge  slowly  if  a  significsuit  high  frequency  compo¬ 
nent  exists.  Furthermore,  when  a  significant  very  low  frequency  component 
(approaching  zero  frequency)  exists,  the  calculated  zero-line  may  still  ignore  an 
excessive  number  of  important  stress  fluctuations. 
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rAnoliteT^ appro achr  ■which  has  gai&ad  iavar  la  aclAvltlas  iacluding 

ocean  wave  hnalysla  and  gust  anaiysis,  involves  defining  a"mlnlmum  included" 
etrain  change.  This  concept  may  be  applied  by  examining  each  crest,  level  Sj, 
together  with  the  trough  preceding  it,  Si.i,  and  the  one  following,  Si+i.  If 
(Si-Si»i)  and  are  both  greater  lhah  the  minimum  to  be  included,  tijen  Sj 

is  defined  to  be  a  valid  crest.  If  that  condition  is  not  met,  then  Si  is  considered  to 
be  a  minor  perturbation  and  is  ignored  in  further  data  processing.  Figure  16  il¬ 
lustrates  the  procedure.  This  definition  of  significant  crests  and  troughs  appears 
to  be  the  best  when  analyzing  records  in  terms  of  alternating  stress  and  mean  stress, 
provided  a  physically  meaningful  truncating  level  may  be  determined. 

In  deriving  a  truncation  level,  It  becomes  necessary  to  examine  a  number  of 
records  in  the  light  of  some  knowledge  of  fatigue.  A  cursory  examination  of  a 
strain  record  may  be  performed  and  the  strain  fluctuation  which  appears  to  be  the 
largest  S^^iax*  rnay  be  selected.  Considerable  fatigue  testing  (Reference  10)  has 
shown  that  fluctuations  smaller  than  will  have  slight  effect  upon  the  fatigue 

life  unless  they  constitute  the  great  majority  of  all  stress  cycles;  examinations  of 
typical  strain  response  caused  by  acoustic  noise  shows  that  fluctuations  smaller  than 
®max/^  constitute  approximately  one-half  of  the  total  number  of  fluctuations.  By 
those  criteria,  then,  a  truncation  level  of  Smax/^  may  be  appropriate. 

However,  more  detailed  study  of  the  typical  strain  records  reveals  that  dis¬ 
carding  the  lower  50%  of  the  strain  fluctuations  invariably  discards  the  mode  of  the 
distribution.  This  is  undesirable  because  it  increases  the  potential  error  incurred 
in  matching  an  empirical  probability  distribution  to  the  measured  one. 

Use  of  a  truncation  level  of  Smax^'^  reduces  the  number  of  discarded  fluctuation 
to  33%  of  the  total,  and  includes  the  mode  of  the  distribution  with  a  small  margin, 

A  larger  margin  may  be  achieved  by  use  of  still  smaller  truncation  levels. 
However,  if  a  level  as  small  as  S^ax/^  is  used,  then  it  is  found  that  a  significant 
number  of  fluctuations  are  dissected  into  a  trivial  alternating  stress  Buperimposed 
upon  a  large  mean  stress,  thereby  jeopardizing  the  alternating  stress  hypothesis. 

The  truncation  level  of  Smax/4  appears  to  be  the  reasonable  compromise,  and 
that  level  has  been  used  throughout  this  report. 


Data  Presentation 


Special  graphic  displays  of  data  are  often  developed  to  facilitate  comparison 
between  theoretical  and  experimental  results.  One  of  the  most  convenient  displays 
of  results  is  one  which  represents  the  results  as  straight  lines  if  the  theory  is  con¬ 
firmed.  Deviations  from  the  straight  lines  are  good  qualitative  measures  of  the 
deviation  of  experimental  results  from  theoretical  results. 

Theory  suggests  a  probability  distribution  of  alternating  stresses  which  ap¬ 
proaches  a  Rayleigh  distribution.  However,  the  truncation  procedure  used  in  the 
data  processing  distorts  the  experimental  results  so  that  they  may  not  be  compared 
directly  to  a  Rayleigh  distribution.  Therefore,  special  graph  paper  was  prepared, 
upon  which  a  truncated  Rayleigh  distribution  plots  as  a  straight  line.  The  truncation 
level  used  in  preparing  the  paper,  0.  9  O'  ,  was  chosen  on  the  basis  that  truncation 
occurs  at  a  level  of  with  determined  from  a  cursory  examination  of 

a  strain  record.  Experience  has  shown  that  the  maximum  determined  from  that  type 
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FIGURE  16  PROCEDURE  FOR  DISCRIMINATING  BETWEEN 
VALID  CRESTS  AND  TROUGHS,  AND  MINOR  EXTREMA 
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of  is  jener ^ly  a  level  which  ia  e«ee<ted  ^  ^^aciately  I  to  2  paaka 

per  thouaanar&ireforerSmax  i®  3-6cr  for  a  Rayietgh  diatributii® 

Smax/4  is  approximately  0.  9  cr 

For  the  inatantaneoua  mean  atreas,  theory  Implies  a  Qaussian  distribution 
The  truncation  procedure  will  exert  some  bias  upon  &e  mean  stress,  since  mo^t 
of  the  stress  fluctuations  which  are  ignored  in  the  data  processing  will  have  small 
mean  stress.  However,  efforts  to  derive  a  procedure  vhich  can  correct  for  thU 
bias  have  not  been  successful.  Because  the  mean  stress  hag  been  assumed  to 
cause  only  second  order  effect  upon  fatigue  life  in  any  event,  it  was  elected  not  to 
pursue  those  efforts  and  to  tentatively  neglect  the  bias.  Therefore,  the  probabilL 
distribution  of  mean  stress  is  plotted  upon  normal  probability  paper. 
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APPENDIX  U1 

TEST  PLAN  FOR  DETERMINING 
EFFECTS  OF  MEAN  STRESS 


Much  of  the  value  of  the  results  of  the  studies  reported  here  depends  upon 
the  verification  of  the  hypothesis  that  a  random  varying  mean  stress,  with  variance 
smaller  than  the  variance  of  the  alternating  stress,  will  have  slight  effect  upon  the 
fatigue  life  of  a  structure.  Therefore,  it  is  considered  appropriate  to  include 
herein  a  possible  test  plan  for  examining  that  hypothesis. 

It  appears  that  the  most  convenient  test  program  for  evaluating  the  relative 
effects  of  mean  stress  and  alternating  stress  will  make  use  of  the  same  natured 
processes  which  tend  to  cause  typical  aircraft  structure  to  be  subjected  to  variable 
mean  stress  and  alternating  stress;  that  is,  response  of  a  several  degree  of  free¬ 
dom  elastic  system  to  broad-band  random  excitation.  A  test  specimen  can  be  de¬ 
signed  to  have  two  well  separated  resonant  frequencies  in  the  frequency  band  of 
interest.  When  both  resonant  modes  are  excited,  the  response  in  the  higher  fre¬ 
quency  mode  will  be  the  alternating  stress  and  the  response  in  the  lower  frequency 
mode  will  be  the  mean  stress.  The  exciter  may  most  conveniently  be  a  random 
vibration  system  with  equipment  for  shaping  the  power  spectrum  of  the  vibratory 
input  so  that  the  relative  response  of  the  two  modes  may  be  adjusted. 

Test  procedure  should  consist  of  first  exciting  the  higher  resonant  mode 
alone.  A  "run"  would  provide  excitation  of  the  higher  frequency  mode  at  a  parti¬ 
cular  mean  square  stress  level  until  failure  is  observed,  thereby  establishing  a 
point  on  a  random  S-N  curve.  Sufficient  runs  should  be  made  to  determine  random 
fatigue  curves  for  each  material  at  negligible  mean  stress.  Succeeding  tests  of 
nominally  identical  specimens  should  then  be  made  with  the  lower  frequency  mode 
excited  simultaneously  with  the  higher  frequency.  The  result  would  be  a  family  of 
curves  of  root  mean  square  alternating  stress  versus  number  of  cycles  of  alter¬ 
nating  stress  to  failure,  for  several  realistic  ratios  of  the  root  mean  square  mean 
stress  to  the  root  mean  square  alternating  stress.  In  order  to  confirm  the  hypo¬ 
thesis  under  trial,  the  differences  between  the  several  curves  of  the  family  should 
be  small;  significant  differences  would  be  cause  for  rejecting  the  hypothesis  and 
showing  instead,  that  the  detailed  form  of  the  stress  power  spectral  density  is 
important  and  must  be  taken  into  account  by  some  method  beyond  the  scope  of  this 
report. 

The  total  test  program  should  include  tests  of  several  materials  and  should 
include  tests  of  both  smooth  and  notched  specimens.  All  specimens  should  be 
instrumented  sufficiently  to  determine  stress  at  the  point  of  failure. 

It  is  urged  that  the  above  test  program  be  initiated  as  soon  as  possible,  so 
that  if  the  results  should  be  favorable  it  will  be  possible  to  take  advantage  cf  the 
other  results  of  this  report. 
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